Dynamical Dirac Mass Generation in the Supersymmetric Nambu— Jona-Lasinio 
Model with the Seesaw Mechanism of Neutrinos 
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The dynamical generation of Dirac mass in the supersymmetric Nambu— Jona-Lasinio ( SNJL ) 
model with the seesaw mechanism of neutrino is investigeted. The right and left handed Majorana 
mass parameters are introduced into the SNJL model; we regard them as external model parameters. 
The question on the origin of these Majorana masses are set aside, and we concentrate on the 
examination of the effect of the Majorana mass parameters on the dynamical generation of Dirac 
mass. The effective potential of the model and the gap equation for the self-consistent determination 
of Dirac mass are derived and solved. We use both the four-dimensional covariant and three- 
dimensional non-covariant cutoff schemes for the regularizations of the effective potential. We find 
there are cases of the first and second order phase transitions with respect to variation of the coupling 
constant of the Nambu— Jona-Lasinio-type four-body interaction of the SNJL model. In the case of 
second-order phase transition, the dynamically generated Dirac mass \4>s\ can arbitrarily be small 
compared with the right-handed Majorana mass parameter \M\ and thus the seesaw condition 
< \4>s\ <C \M\ can be satisfied by a fine tuning of the coupling constant, while at the first-order 
case it seems very difficult and/or "unnatural" to satisfy the condition. The numerical results do 
not depend on the difference of the cutoff schemes qualitatively. 

PACS numbers: 11.30.Pb,12.60. Jv,12.60.Rc,14.60.St 

Recent experimental observations of flavor oscillations confirmed that neutrinos have small masses ( the mass of 
tau neutrino, < 18 McV ) compared with other particles inside their same generations. This fact motivates us to 
consider the seesaw mechanism of neutrinos [1-3] more seriously, because the mechanism seems to provide the unique 
explanation on tiny masses of neutrinos. The seesaw mechanism states that Dirac masses ( for example, ~ 300 GeV ) 
of neutrinos will be suppressed by very heavy Majorana masses ( for example, 2 x 10 16 GeV ). The electroweak gauge 
theory of the standard model is the most important success of modern particle physics, which gives very accurate 
agreements with various experimental results, and still it will give us the horizon of particle phenomenology In 
the electroweak symmetry breaking, the top condensation model [4,5] obtains much attractions from us until now, 
because it is one of possible candidates toward beyond the standard model. The top condensation scenario uses a 
generalization of the Nambu— Jona-Lasinio ( NJL ) model [6] to include the electroweak gauge symmetry, and it 
considers a fermion-antifermion condensate. Recently, Antusch et al. gave a theory of an NJL-type four-body contact 
interaction model with introducing a right-handed Majorana mass parameter [7], In that work, the authors derived a 
gap equation for the dynamical generation of a Dirac mass term of neutrino, and solve them numerically. Moreover, 
they argued the standard model with right-handed neutrinos as a low-energy effective theory of their model under 
its condensation scale A, and performed a renormalization group analysis of model parameters of the renormalizable 
effective theory. It is an interesting problem for us to examine a supersymmetric version of a four-fermion model 
combined with a right-handed Majorana mass term for neutrino. This is the motivation of this paper. We employ the 
supersymmetric Nambu— Jona-Lasinio ( SNJL ) model as our model Lagrangian. Historically, the SNJL model was 
introduced for the investigation of dynamical chiral symmetry breaking in supersymmetric quantum field theory [8] , 
and it is extended to construct a supersymmetric version of the top condensation model [9] . The SNJL model was 
also used to investigate the SUSY BCS-type superconductivity [10,11]. 



We introduce the following one-flavor SNJL model with including the Majorana mass parameters: 
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Here, <&+ and <&_ are right and left handed chiral superfields of neutrino, respectively. We take the convention of 
spinor algebra by the book of Wess and Bagger [12]. M and m are right- and left-handed Majorana mass parameters 
( in general, given as complex numbers ), respectively. The term of the coupling constant G in (2) will be used as a 
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four-body contact interaction similar to the case of the ordinary Nambu— Jona-Lasinio model [6]. In four-dimensional 
spacetime, G has mass dimension [mass] -2 and the theory is non-renormalizable. In this paper, we perform the 
four-dimensional covariant and three-dimensional non-covariant cutoff regularizations, and a choice of a numerical 
value of the cutoff will be regarded as an introduction of the limit of the applicability of the theory by our hands. The 
electroweak gauge symmetry is not considered in this model for the sake of simplicity, hence it is a drastically simplified 
model prepared for our examination on the dynamical generation of a Dirac mass under including the Majorana mass 
parameters. In the seesaw mechanism of neutrino masses, a right-handed Majorana mass of neutrino should take a 
huge value compared with its Dirac mass. The mechanism of the generation of Majorana mass might be given in the 
different framework from the dynamical generation of Dirac mass ( for example, by a vacuum expectation value of a 
Higgs field ), hence we regard M ( and also m ) as a model parameter. Of course, a simultaneous consideration of the 
dynamical generations of both the Dirac and the Majorana mass terms can be possible, by using the similar framework 
of the supersymmetric theory of superconductivity [10,11]. Here, we assume the generation of the Majorana mass 
terms happens outside, a possible underlying theory, of our model. In this paper, we just focus on the examination 
on some effects of the existence of the Majorana mass parameters in the dynamical generation of a Dirac mass term. 
Our attitude to the problem is similar to the work of Ref. [7], and we can say that our theory is a SUSY version 
of Ref. [7] , though we do not consider the electroweak symmetry breaking. Usually, it is believed that the existence 
of the SUSY is favorable to solve the problem of gauge hierarchy. Because our model is determined as a one-flavor 
model, we do not consider the Kobayashi-Maskawa-like mixing matrix in our theory. By using this model, we derive 
the effective potential and the gap equation for a self-consistent Dirac mass. 

By using the method of SUSY auxiliary fields of composites [8-11], the model Lagrangian will be converted into 
the following form: 



C = 
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Here, S and S' are Lagrange multiplier multiplets to keep the following definition of the composite collective chiral and 
antichiral superfields H = G$+$_, = G$>' + &_, respectively. Due to the non-renormalization theorem of super- 
symmetry, the ordinary SNJL model with keeping Af = 1 SUSY cannot break the chiral symmetry dynamically [8,9]. 
We consider it is also the case in our model. In order to generate a Dirac mass term in our theory, we have introduced 
a SUSY breaking mass A in (5). Expanding C in terms of component fields, eliminating the auxiliary fields of chiral 
multiplets through their Euler-Lagrange equations, and assuming that the collective fields are spacetime- independent, 
one finds the Lagrangian in the following form: 



C = - 



\<t>sf 



G 

-(|<fe| 2 + 4|M| 2 + A 2 )<^+ - (|0s| 2 + 4|m| 2 + A 2 )0U 
+{2M(j)s + 2m t 4)</>+</'- + (2Mt(4 + 2m<j)s)^_4> + 



(6) 



At the case |M| ^ \m\, the Lagrangian has a skew structure in the right-left symmetry space, and thus the contributions 
of the right and left handed particles in quantum corrections of a loop expansion will have differences. The generating 
functional is given as 



V^ + V^ + V^_Vi}_V4> + V(j)\V<j)_V<j)lV(j)sV(j)\ 
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where, the matrices fig. and Sip are given as follows: 



() = , fc 2 -|0 s | 2 -4|M| 2 -A 2 2{M<f>s + m^ s ) , , . 

iLB - \ 2(Mt4 + m0s) fc 2 -|0 s | 2 -4|m| 2 -A 2 ' ' 



( 7 5 = 7°7 1 7 2 7 3 ) while the definitions of several fields in our Lagrangian are 

* = (iPmr,*Pml) T , ^mr = (■0+,'0+) T , iPml = (ip-,$-) T , 

* = {lf>MR,1pML), tpMR = (-i> + ,-i> + ), IpML = {-l)) -,-$-) , II = ((f) + , 0_ ) T . (10) 

Here, i/'Mit and V'ml are right- and left- handed Majorana fields, respectively. T denotes transposition. ^ can be 
called as an eight-component Nambu-notation-field [13] defined in terms of the Majorana fields. ( Of course, it is 
equivalent to the formalism of eight-component Nambu-notation given by Dirac fields which was used in relativistic 
theory of superconductivity [14,15]. ) Because we have three complex mass parameters 

4>s = \4>sW 9s , M=\M\e l8M , m=\m\e l0m , (11) 

we can absorb only two of these phases 0g, 9m and 6 m by a redefinition of fields ip+ and ip_, and one phase degree 
of freedom remains. 

We employ the steepest descent approximation for the integration of the collective fields <fis and <fr s . Therefore, we 
obtain the effective action in the following form: 

T e// = J dAx (-^-) +2ilnDctf2 B -ilnDetfi F . (12) 
The effective potential will be obtained after the diagonalizations of the matrices f2 b and as 

V l^l 2 f d ' k i {^-E B {k)){k + E B {k)){k^E B {k)){k + E B {k)) 

ft G V (2tt) 4 11 (k a -El(k))(k + El{k))(k -E F {k))(k a + EZ{k)Y { ' 

( the diagonalization of f2_p has been done by the same method given in Ref. [14] ) where, the " quasiparticle" excitation 
energy spectra become 



E B (k) = Jk 2 + \cj) S \ 2 + 2\M\ 2 + 2\m\ 2 + A 2 T 2\/(|M| 2 - |m| 2 ) 2 + |(fe| 2 (|M| 2 + |m| 2 + 2\M\\m\ cosO), (14) 



E£(k) = \Jk 2 + \4> s \ 2 + 2|M| 2 + 2|m| 2 t 2^/(\M\ 2 - |m| 2 ) 2 + |<fe| 2 (|M| 2 + |m| 2 + 2\M\\m\ cosO), (15) 

e = 20s + e M + e m . (16) 

The masses appear in these spectra show complicated structures, though the Lorentz symmetry is still kept in our 
theory. ( These spectra have some similarities with the mass spectra of top/bottom-stop/sbottom appear in the 
Minimal Supersymmetric Standard Model (MSSM) [16-18]. ) We show several limiting cases of these spectra. Under 
| M | = |m|, these spectra become like 

E B (k) = v /fe 2 + (|0 s | T 2|M|) 2 + A 2 , (17) 
££(*) = ^fe 2 + (|0 s |T2|M|) 2 (18) 
at 6 = 0. On the other hand, our E± and E± become 



E B (k) = Jk 2 + \<f> s \ 2 + 4\M\ 2 + A 2 , (19) 



E£(k) = V /fc 2 + |0 s | 2 + 4|M| 2 (20) 
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at O = 7r and they coincide with the dispersion relations of the zero-chemical-potential case in the SUSY BCS 
superconductivity [10,11]. In fact, the BCS-type superconductivity is a dynamical generation of a specific choice of 
the right- and left- handed Majorana masses. At \M\ ^ 0, \<j>s\ ^ 0, |m| = 0, 



El 



(k) = ^k 2 + {\M\ T V|M| 2 + |0 S | 2 ) 2 ^A2, 



E£{k) = ^k 2 + (\M\ T V\M\ 2 + \ct>s\ 2 ) 2 , 
and thus, in the case \M\ » \<j)s\ > with \m\ = ( the seesaw condition ), the spectra become 



Ef(k) 



2 , \M 



41 M 



+ A 2 , 



E B (k) = ^/fc 2 + 2|(/) S | 2 +4|M| 2 + A 2 , 
E F + {k) = Jl 
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£ F (fc) = ^/fc 2 + 2| ( /) S | 2 +4|M| 2 . 



(21) 
(22) 

(23) 
(24) 
(25) 
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Hence, the seesaw mechanism suppresses the mass of the branch of E+(k) just the same with the ordinary seesaw 
theory [1-3]. In our theory, the seesaw mechanism occurs also in the scalar sector. The choice |m| = reduces the 
number of mass parameters in our theory, and thus the phase degree 6 disappears from our theory. 

After performing the Wick rotation in the momentum integration J d A k of V e ff and introducing the covariant 
four-momentum cutoff, the effective potential is found to be 



V 



1 



eff 



\<i>s\ 2 



G 16tt 2 



2A 2 A 2 + A4ln( 1 + a +/ A ^ 1 + a -/ A2 ) 



Here, A denotes the four-momentum cutoff. a± and (3± correspond to the squares of the masses appear in E±, E±. 



(l + /3+/A 2 )(l + /3_/A 2 ) 
-a\ ln(l + A 2 /a+) - a 2 _ ln(l + A 2 /a_) + /? 2 ln(l + A 2 / f3+) + /? 2 ln(l + A 2 //3_) 

a± = (Ei(k = 0)) 2 , 0± = {Ei{k = 0)) 2 . 

tes the four-mor 
At A 2 3> a± , fi± , one obtains 

v l4)sl \ 1 

Veff = — + 

We also examine our V e f / in the three-dimensional non-covariant cutoff scheme: 

Veff = 



(27) 
.., E£. 



A 2 A 2 A 2 A 2 

4A 2 A 2 - a 2 . In — - a 2 _ In — + /3 2 In tL. + f_ m ±- 
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(28) 
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-a + ln ■ 
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v/A 2 +/3_+A' 



G 8vr 2 

-2A(A 2 + 13+ f' 2 + A/3 +v /A 2 + /3+- 2A(A 2 + f3- f 2 + A/3_ ^/WTp 
VA 2 + a+ + A 



a 



• (29) 



The energy spectra of particles we have obtained are Lorentz symmetric, while the three-dimensional cutoff scheme ( 
usually be employed in a finite-temperature Matsubara formalism ) breaks the Lorentz symmetry explicitly. Therefore, 
the self-consistent gap equation for |</>g| is found from the stationary condition of the effective potential 



dv a 



be 



d\4>s\ 



to 
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G 
1 



d\<p 



A 2 
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(30) 
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in the covariant cutoff scheme, while 



8V, 



eff 



d\<t>s 



2|fo| 1 
G 4tt 2 

da 



<9aq 



^ AVA 2 + a+ - AVA 2 + [3+-a+ In — ^ + (3+ In ^= ^ 



f A VA^T - Av^TTT - a In V^L±^ + ^ ln VEIj£±A\ 



in the three-dimensional cutoff regularization. Here, the derivatives appear in the gap equations become 



dot± 



d\<j>s\ d\<j>s\ 



= 2 



\<t>s\ ( 



IT 



|M| 2 + |m| 2 + 2|M||m|cose 



^(|M| 2 - H 2 ) 2 + \4> S \ 2 (\M\ 2 + H 2 + 2\M\\m\ cosG) 



(31) 



(32) 



Hereafter, we examine the case \m\ = and thus the phase <d disappears from V e ff and the gap equations. We 
used the numerical package Mathematica ver. 6 for our calculation. At the case \m\ = 0, the determination equation 
for the critical coupling in the covariant cutoff scheme is obtained as follows: 



G c 



4tt 2 



(A 2 + 4|M| 2 ) ln(l + - 4|M | 2 ln(l + ^) ' 



(33) 



From this results, we obtain lim| M |^ G cr = 4tt 2 /{A 2 ln(l + coinsides with that of the ordinary SNJL model [8]. 

Figure 1 shows the critical coupling G cr as a function of cutoff A under the unit A = 1. G cr shows the well-known 
dependence on A, and becomes small at a large A. G cr depends on \M\ quite sensitively, and when \M\/A becomes 
large, the critical coupling also becomes large and the dynamical generation of a finite VEV of \4>s\ will be suppressed. 
This stems from the fact that the denominator of Eq. (33) approaches to zero under the limit \M\/A — ► +oo. The 
critical coupling in the three-dimensional cutoff scheme is found to be 



G c 



2n A 



A^A 2 + 4|M| 2 - AVA 2 + 4|Af| 2 + A 2 



-(4|M| 2 + A 2 )ln 



v /A 2 + 4|Af| 2 + A 2 +A 
V4IMI 2 + A 2= 



4|M| 2 ln 



y/A 2 +4\M\< 

2\M\ 



A 



(34) 



This expression of G cr of the three-dimensional non-covariant cutoff with \M\ — coinsides with the results of 
Refs. [10,11] of the zero-density case. Figure 2 shows G cr of the non-covariant cutoff scheme as a function of A. The 
result of Fig. 2 is qualitatively the same with that of Fig. 1. In the covariant cutoff scheme, the second derivative of 
V e ff with respect to \(f>s\ at the origin will be given by G cr of Eq. (33) as follows: 



8 2 V f 



eff 



\<Ps\=0 



\G G cr J 



(35) 



The numerical results of V e // of Eq. (27) by the covariant cutoff scheme are shown in Figs. 3, 4, 5 and 6, while that 
of Eq. (29) of the non-covariant cutoff scheme are given in Figs. 7, 8 and 9. All of our numerical calculations of V e // 
are performed under a fixed SUSY breaking scale A: A/A is fixed to 100, and we use the energy unit A = 1. The 
global minimum of V e ff is quite sensitive to both G and \M\, and thus we regard them as the relevant parameters 
to control a numerical value of self-consistent mass \<f>s\- A too large G gives an unphysical solution like \(ps\ > A. 
Hence we will choose a numerical value of G close to G cr . In our numerical calculations, we find there are cases of 
first and second order phase transitions they take place under the variation with respect to G. In the case of first 
order transition shown in Figs. 3, 4, 5, 7 and 8, they are depicted by several numerical examples of G, the global 
minimum of V e f f usually locates at \<$>s\ ~ C(| Af |) and such a \<$>s\ 1S too large for satisfying the seesaw condition 
< 1 4>s | «C \M\ in our model. In this case, when we take a small value of G to try to obtain a \<j>s\ "C \M\, V e ff 
suddenly changes the global minimum from \<f>s\ ^ into its origin \(j>s\ = 0. Hence it is difficult to obtain a small and 
finite VEV of \</>s\ which satisfies the seesaw condition. When \M\ is large, G cr also becomes large, and the energy 
scale of the tree level part \<fis\ 2 /G becomes large for obtaining a finite VEV of \4>s\- See the differences of energy 
scales of Figs. 3, 4, 5 and 6. The reason of \cps\ ~ C(|A4"|) stems from this fact under the first-order transition. 
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On the other hand, in the case of second order transitions given in Figs. 6 and 9, a global minimum of V e ff can 
take an arbitrary small value of \(f>s\ compared with the right-handed Majorana mass \M\ by choosing an appropriate 
numerical value of G, and thus the seesaw condition < \(f>s\ <C \M\ can be satisfied. In this case, the theory 
approaches to a critical point with a continuously vanishing \<f>s\- 

Our model becomes the ordinary SNJL model which has no Majorana mass term, at the limit \M\ — > with fixed 
A and A, while the non-SUSY NJL model will be obtained by A/A — > oo. The ordinary SNJL shows second-order 
transition under the variation of G. Therefore, there is a boundary of the first and the second order phase transitions in 
the parameter space of the ratio \M\/A. For example, we found in our numerical calculation, the boundary between 
the second and the first order transitions locates at the region 0.5 < \M\/A < 0.55 under the parameter choice 
\m\ = 0, A = 1, and A = 100 in the covariant cutoff scheme: In this case, the phase transition becomes second-order 
at |M|/A < 0.5, while it becomes first order at |M|/A > 0.55. ( Similarly, we found the boundary in the region 
0.45 < \M\/A < 0.55 under the choice |m| = 0, A = 1, and A = 100 in the case of non-covariant cutoff scheme. ) 
The boundary could be understood as 



Therefore, we obtain the conclusion that the right-handed Majorana mass \M\ should take a value much smaller than 
the SUSY breaking mass scale A ( in the case A 3> A ), which is considered to have a TeV energy scale mass in 
modern particle phenomenology, to satisfy the seesaw condition \M\ » \<j)g\ > in our model (1). In the case of 
second-order transition, we can choose G close to G cr ( a fine-tuning procedure [5] ), while in the first-order case it 
is impossible to do the procedure to obtain a finite VEV of \<f>s\ which satisfies the seesaw condition. At the origion 

of the effective potential, the second derivative g i ^jl of (35) is always negative under the second-order transition, 
while it will take positive values under the first-order transition. This means G > G cr is always satisfied for obtaining 
a finite VEV of \<fts\ m the case of second-order transition, while V e ff can have a finite VEV \(f>s\ ^ even if G < G cr 
under the first-order transition. ( Therefore, G cr lost the precise meaning under the first-order transition case. ) 
These results of V e f / we have obtained were also confirmed in our numerical evaluation of the gap equations (30) and 

(31), and also by the behavior of second derivative §rafrj at the origin \cps\ = of the cases of Figs. 3-9. The author 
also examined the A-dependence of V e ff, by choosing A/A — 500. Still, the case |M|/A = 0.1 shows the second-order 
transition, while the case \M |/A = 1 gives the first-order transition. Hence the order of transition is determined by 
the relation between \M\ and A under \M\, A <C A, as indicated in (36). Figure 10 shows the first-order transition of 
the example A = 1, |M| = 1000 and A = 1000000. 

Our conclusions are summarized as follows: 

• We have shown that the dynamical Dirac mass with the seesaw mechanism of neutrino can be described by our 
simplified model. 

• When \M\ > A or \M\ ~ A with A ^ A, a self-consistent solution of \<fis\ will take a numerical value of the 
order 0(|M|). In this case, it is difficult to satisfy the seesaw condition < \<j>s \ "C \M\. 

• When \M\ is much smaller than A with A » A, we can obtain an arbitrarily small \<ps\ by variation with respect 
to G. In this case, it is possible to satisfy the seesaw condition. 

The final comment. Recently, the Nambu— Jona-Lasinio and Gross-Neveu models have been derived from string 
theory with intersecting D-branes [19-21]. It is intersting for us whether the Majorana mass parameters we have 
considered here can be derived from such an attempt, and they ultimately have their origins in string theory, or 
not. In another point of view, effective actions obtained from compactifications of supcrstring theory has a common 
feature, no matter what method of compactification one uses [22,23,24]. While, the MSSM can be regarded as the 
low-energy effective theory of the top-condensation model. Hence, our ultimate goal is summarized in the following 
scheme ( see, also [25] ): 



2|M|/A - 1. 



(36) 



Supcrstrings -> SNJL + Majorana -> MSSM. 



(37) 



* Present address: Department of Physics, University of Texas at Austin,tadafumi@physics. utexas.edu 
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FIG. 1: The critical coupling constant G cr of the covariant cutoff scheme as a function of cutoff A. We choose A = 1. 

FIG. 2: The critical coupling constant G cr of the non-covariant cutoff scheme as a function of cutoff A. We choose A = 1. 

FIG. 3: The effective potential V e // ( obtained by the covariant cutoff scheme ) of the first-order phase transition as a function 
of the dynamical Dirac mass parameter \<j>s\- We set the parameters as A = 1, \M\ = 10, A = 100, \m\ = 0. In this case, we 
obtain G cr A 2 = 17.1989 from Eq. (33). The second derivative ( i.e., the curvature ) of V e // at the origin is negative at the 
case GA 2 = 17.2, while it is positive at the cases GA 2 = 16.1, 16.0, 15.9. 

FIG. 4: The effective potential ( obtained by the covariant cutoff scheme ) of the first-order phase transition as a function of 
\<f>s | • The parameters A, |M|, A and \m\ are set as the same with Fig. 3. The positiveness of curvature of V e // at the origin is 
clear under the energy scale of this figure. 

FIG. 5: The effective potential ( obtained by the covariant cutoff scheme ) of the first-order phase transition as a function of 
\4>s\- We set the parameters as A = 1, \M\ = 1, A = 100, \m\ = 0. 

FIG. 6: The effective potential ( obtained by the covariant cutoff scheme ) of the second-order phase transition as a function 
of \4> s \. We set the parameters as A = 1, |M| = 0.1, A = 100, \m\ = 0. 

FIG. 7: The effective potential ( obtained by the non-covariant cutoff scheme ) of the first-order phase transition as a function 
of \4> s \. We set the parameters as A = 1, \M\ = 10, A = 100, \m\ = 0. 

FIG. 8: The effective potential ( obtained by the non-covariant cutoff scheme ) of the first-order phase transition as a function 
of \4> s \- We set the parameters as A = 1, \M\ = 1, A = 100, \m\ = 0. 

FIG. 9: The effective potential ( obtained by the non-covariant cutoff scheme ) of the second-order phase transition as a function 
of \4> s \. We set the parameters as A = 1, \M\ = 0.1, A = 100, \m\ = 0. 

FIG. 10: The effective potential ( obtained by the covariant cutoff scheme ) of the first-order phase transition as a function of 
\(f>s\- We set the parameters as A = 1, |M| = 1000, A = 1000000 and \m\ = 0. 
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